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We show a model of the slow expansion, in which the scale invariant spectrum of curvature 
perturbation is adiabatically induced by its increasing mode, by applying a generalized Galileon 
field. In this model, initially e <C — 1, which then is rapidly increasing, during this period the 
universe is slowly expanding. There is not the ghost instability, the perturbation theory is healthy. 
When e ~ — 1, the slow expansion phase ends, and the available energy of field can be released and 
the universe reheats. This scenario might be a viable design of the early universe. 
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I. INTRODUCTION 



The observations imply that the primordial curvature 
perturbation is scale invariant. Thus how generating it 
has been still a significant issue, especially for single field. 
The curvature perturbation on large scale consists of a 
constant mode and a mode dependent of time [l| . When 
one of which is dominated and scale invariant, the spec- 
trum of curvature perturbation will be scale invariant. 
When the scale factor is rapidly changed while e is nearly 
constant, the constant mode is responsible for that of in- 
flation while the increasing mode is for the 
contraction with matter [(| , Q , [H , both are dual Q. 

In principle, the increasing mode of metric perturba- 
tion, which is scale invariant for e > 1 |) or e < -1 
[lfjj . might dominate the curvature perturbation. The 
constant mode of metric perturbation is same with the 
constant mode of curvature perturbation. The duality of 
scale invariant spectrum of metric perturbation has been 
discussed in [ll[ , [l^l , [H| • The evolution with e >• 1 is 
the slowly contracting, which is that of ekpyrotic uni- 
verse [Hj]. While e -C — 1 gives the slow expansion [To| . 
which has been applied for island universe [HI- In cer- 
tain sense, in Ref.pjJ it was for the first time observed 
that the slow expansion might adiabatically generate the 
scale invariant spectrum of curvature perturbation, see 
[l6j for that induced by the entropy perturbation. 

When the available energy of field is released, the slow 
expansion phase ends and the universe reheats. Thus the 
slow expansion might be a viable scenario of the early 
universe. In principle, when e is constant, whether the 
increasing mode of the metric perturbation can be inher- 
ited by the curvature perturbation depends of the physics 



that the evolution of e < emerges only for a period, the 
phantom field might be only a simulation of a full theory 
without the ghost below certain physical cutoff [22j . 

Recently, the cosmological application of Galileon, 
[13|,[13|> or nontrivial generalization |25 | . l26l. l27l. has 
acquired increasing attentions [28j . [29| . [30| , [31| , |32| . It 
has been found for generalized Galileon that e < can 
be implemented stably, there is not the ghost instability. 
We, in this paper, will show a model of the slow ex- 
pansion given in [19( , by applying a generalized Galileon 
field. In this model, the perturbation theory is healthy, 
the scale invariant curvature perturbation is given by it- 
self increasing mode, which can be consistent with the 
observations. As will be argued, this in certain sense 
validates the argument and calculations in [Toj.[l9j 

The models of early universe, builded by appl yin g 
generalized Galileon, have been studied. In Ref . [26| . 
the inflation model is implemented by using generalized 
Galileon field. However, here what we discuss is an al- 
ternative to inflation. There is a slightly similar scenario 
in [33j . However, in [33| . the adiabatic perturbation is 
not scale invariant, the scale invariant curvature pertur- 
bation is obtained by the conversion of the perturbations 
of other light scalar fields. Here, we will see how the 
adiabatic perturbation is naturally scale invariant. 



II. AS A GENERAL RESULT 

We begin with a brief review on the slowly evolving 
model in [jjjj. The quadratic action of the curvature 
perturbation 1Z is 



ited by the curvature perturbation depends ot the physics r a ^Q , /2 

around the exiting 17|. However, when e is rapidly $2 ~ / drid x — — \1Z — c s (dlZ)' 



changed, the thing is altered, see [18[ for that of the slow 
contraction. During the slow expansion, the scale invari- 
ant curvature perturbation can be naturally induced by 
its increasing mode 19], or its constant mode [20j.(2lj. 

The perturbation mode can leave the Hubble horizon 
during the slow expansion requires e < [HI, [HI, or a pe- 
riod after it is required to extend the perturbation mode 
out of the Hubble horizon [20]. Thus in [l0|,[l9j], the 
phantom was applied for a phenomcnological studying. 
However, there is a ghost instability. Thus it was argued 



(1) 



which is actually general for single field, like P(X, ip) [34j . 
generalized Galileon [25[ , [2(| , [3{| , and the modified grav- 
ity [lH , [37} ■ Q and c 2 are generally different for different 
models. However, Q > and c 2 s > should be satisfied 
to avoid the ghost and grad ient instabilities. 
The equation of K is [11 ,[11 



+ elk 2 Ufc = 0, 



(2) 



2 



after defining Uk = 
7], z = a^/2MpQ/c s 
constant c 2 . 
by Eq.© is 



zTZk, where ' is the derivative for 
. We here only care the case with 
When fc 2 <C z" /z, the solution of 1Z given 



1Z ~ C is constant mode (3) 

— is changed mode, (4) 

z z 

where D mode is increasing or decaying dependent of 
different evolutions. 

The scale invariance of 1Z requires ~ 2 „^ 2 , which 
implies 



1 



for constant mode 



(5) 



V* - V 

(r)* — rj) 2 for increasing mode (6) 



has to be satisfied, where initially rj -C — 1. In certain 
sense, both evolutions are dual @ . The results will be dif- 
ferent if c 2 s is changed, however, which we will not involve 
here. In principle, both a and Q can be changed, and to- 
gether contribute the change of z. However, only one 
among them is changed while another is hardly changed 
might be interesting, e.g. the inflation, given by (5), or 
the contraction dominated by the matter, given by ((6|), 
in which a is rapidly changed while Q is hardly changed. 

However, the case can also be inverse. When Q is 
rapidly changed while a is hardly changed, the scale in- 
variant spectrum of curvature perturbation can also be 
induced by cither its constant mode [l8[ , [20| , [2l| , given 
by (5), or its increasing mode [Hj], given by Though 
both cases give the scale invariant spectrum, both pic- 
tures are distinct. In general, for the picture in [19j, 
initially |e| 3> 1, which then is rapidly decreasing, the 
slow evolution of the scale factor ends when |e| ~ 1. 
While for that in QJ] , [H , HH , initially |e| < 1 which 



then is rapidly increasing. In addition, for p| , [20]| , [2lj , 
during the slow evolution, the perturbation mode is actu- 
ally still inside the Hubble horizon. Thus a period after 
it is required to extend the perturbation mode out of 
the Hubble horizon, while in |19j | , the perturbation mode 
can naturally leave the Hubble horizon during the slow 
evolution. There is also not the problem pointed in [40| . 

Here, we will discuss that in [Ijj]. We have generally 
Q = e for single field action P(X,<p) [H|. While the 
case is slightly complex for generalized Galileon [25| , [26| . 
However, as will be showed in following section, we actu- 
ally have Q ~ |e|. 

Thus Q = |e| will be set for general discussions in the 
following. In principle, |e| is dependent of a. However, 
it can be observed that a is nearly constant for |e| 3> 1. 
Thus for ([5]), we have 



\e\~Ai(U-t) 4 



(7) 



since rj ~ t, where A* is dimension. The Hubble 
parameter is given by 



H 



1 



A*(t,-i) 6 



(8) 
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FIG. 1: The evolutions of a, the Hubble horizon and the 1Z 
horizon during the slow expansion given by Eq.([7|. a» = 10 
is set. During this phase, due to the rapidly change of H and 
H f ree ze, the perturbation mode initially inside both horizons, 
i.e. A ~ a <C 1/H freeze *C 1/H will naturally leave the 1Z 
horizon, i.e. A ~ a > 1/ H freeze, and then the Hubble horizon, 
i.c.A ~ a > 1/H. 



Thus a is given by 



ln[^ 



1 



Ai(U-t)+ |e| 



(9) 



When initially A*(Tj* — t) 3> 1, i.e.|e| 3> 1, the evolution 
corresponds to the slow expansion for e <C — 1, or the 
slowly contraction for e > 1, since a/a* ~ 1. The slow 
evolution ends when A*(i* — t) ~ 1, at which |e| ~ 1. 



When k 2 



' /z, the perturbation mode is leaving the 



horizon, and hereafter it freezes out. This horizon might 
be called as the 7Z horizon 



1/Hfr 



v* - V- 



(10) 



Thus the physical 7Z horizon is a/H freeze — t* — t, While 
the Hubble horizon is 1/H given by Eq.(j8]). Here, the 
evolutions of the 1Z horizon and the Hubble horizon are 
different. While when a is rapidly changed and |e| is un- 
changed, e.g. inflation, both evolutions are almostly same. 
The reason is that for inflation, z" / z ~ a" /a, thus 



1/Hfr 



z 
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a" 



i/n, (ii) 



while here a is constant and |e| is rapidly changed, we 
have not z" / z ~ a" /a. 

When fc 2 ^> z"/z, i.e. the perturbation is deep inside 
the 1Z horizon, Uk oscillates with a constant amplitude. 
The quantization of Uk is well defined for Q ~ |e| > 0, 
which gives its initial value. The evolutions of a, 1/H and 
a/H freeze are plotted in Fig.l for the slow expansion. It 
can be found that the perturbation mode firstly leaves 
the 1Z horizon, after which it is freezed out, but it is still 
inside the Hubble horizon. However, since the Hubble 
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horizon is decreasing, after a while the perturbation mode 
will be inevitably extended outside it, and become the 
primordial perturbation on super Hubble scale. 

When k 2 -C z"/z, the amplitude of perturbation spec- 



trum is V 



1/2 
K 



Thus 



-H 2 



(12) 



determined by 



^/M + ^_ X 2 + ^ H ipx) ip 



+ 3 



M'< 



4 V /M 



M 



M 1 



M"' 



where Q ~ |e| is applied. The perturbation is given by 
the increasing mode Q , because a is hardly changed and 
|e| is decreasing. When |e| ~ 1, the change of a begins to 
become not negligible. Though |e| is still decreasing, a 
is increased exponentially. Thus this mode will become 
the decaying mode at certain time tf ~ 0(i*) shortly 
after |e| ~ 1. In principle, the spectrum of 1Z should be 
calculated around tf. Thus 



1/2 

n 



Mi 



(13) 



The universe reheats around or after tf, and hereafter the 
perturbation is dominated by its constant mode, until it 
enters into the Hubble horizon during the radiation or 
matter domination. |e/| ~ 1 brings A^(i* — tf) 4 ~ 1. 
Thus Eq.(fT5]l becomes 



,1/2 

n 



A, 



Mr 



(14) 



which is general result of the slow evolution in i.e. 
the evolution of |e| follows Eq.© and c 2 is constant. 



III. A GALILEON DESIGN OF SLOW 
EXPANSION 

Here, we will detailed show a model of the slow ex- 
pansion given in (To|.fl9ll . While the scenario of the slow 
contraction given in [19| is slightly alike with that in [l8[ , 
which might be studied in detail elsewhere. 



A. The background 

We consider a generalized Galileon as 

L^-e 4 ^X + -^X 3 -^X 2 B^ (15) 



M< 



where A4 is the energy scale. Here, the sign before 
£ 4<p/m x i s negative. However, as will be showed that 
this model has not the ghost and gradient instabilities, 
since Q > and c 2 > 0. The evolution of background is 



m 2 M 2 = - e 4v/M X 



5 



X 6 



6 



Ml X^H. (17) 

We require that initially e <C — 1, and behaviors as 
Eq.0. This can be found by requiring e 4v ^ M X 
in Eq.jnj}. This gives 
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1/4 



1 



m{u -ty 



Thus 



Thus 



M 



(u-ty 



H 



M 7 M 2 M 2 {U 



ty 



(18) 



(19) 



(20) 



is induced. Thus for MMp ~ A%, Eq.® is obtained. 
This gives Eq.([7]), which is just required evolution. 

Eqs. flTB)) and (fTT)) are numerically solved in Fig. 2 and 
Fig. 3. We can see that Eqs. (Tl9t and ([20)) can be highly 
consistent with accurate solutions for a long range of 
time. The significant deviation only occurs around tf ~ 
We might think that the slow expanding phase 
ends when the significant deviation appears, and the re- 
heating begins. However, it might be possible that the 
reheating of universe begins some time after the signifi- 
cant deviation occurs, since the perturbation generated 
during this period only are the perturbation on small 
scale, which has not to be scale invariant. 

Eqs.fUS]) and ® implies HtpM < X, Hip/M 3 < 
e 2tp l M , and Hip < ip, since 



H 



1 



1 



(t* - tf (U - t) 



(21) 



for |e| > 1, \.e.^MM P {U - t) > 1. Thus Eq.® is 
approximately 



= 4 V /A4 



(22) 



for \fMMp~{t* - t) > 1. It can be found that Eq. ([22j) is 
consistent with Eqs. p^|) and (jT^J) - Thus the equation of 
the perturbation Sip of ip is 
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FIG. 2: The evolution of ip with respect to the time. The 
initial values of p and ip are required to satisfy Eqs. (|18[) and 
((19) , respectively. The parameter M = O.OlA/p. The dashed 
line is that of Eq. (|19[l . The inset is that around tf ~ 0(t*). 
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FIG. 4: The evolutions of the amplitude of curvature per- 
turbation for different k with respect to the time. The green 
and black lines are that with different k. Here, the time axis 
is rescale as Mt for the convenience of numerical calculation, 
t is that in Fig.2 and Fig.3, M = 0.01. 
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FIG. 3: The evolutions of a and H with respect to the time. 
The red line is that of H. The black line is that of a, while the 
black dashed line is that of Eq.((9). The inset is that around 
t f ~ 0(U). 



When Eqs. (fT8|) and (|19j) are considered, the solution is 

Sip ~ (t* — t) , is decaying mode (23) 
or l/(t* — t), is increasing mode. (24) 

The decaying mode is negligible. The increasing mode 
is dominated. Thus Sip ~ -jgj. Thus for MAt » 1, 
Sip <C Aip. Thus if initially Sip -C ip is satisfied, it will 
be valid all along. When the time arrives around tf, 
Eq.(|2"T]) will be not right. Thus Eq.([2"2|) can not be found. 
This explains why there will be significant deviation for 



Eq.(|l9]l around tf. 

There might be other fluids, However, their energies 
generally do not increase, since the expansion is slow. 
Thus for |e| > 1, i.e.\/MM P (U - t) > 1, the evolution 
of background, given by Eqs. ([TT)f and (|2U| . is stable. 



B. The curvature perturbation 



1Z satisfies Eq. 
culations of Refs.[2C 



We follow the definitions and cal- 
HBt Here, the generalized Galileon 



FIG. 5: The spectrum of curvature perturbation at differ- 
ent times with respect to k. The black dashed line is initial 
spectrum. The short dashed, long dashed and solid orange 
lines are the spectra at different times, respectively. There is 
a cutoff k cu t ff ~ 5 x 10 , below which the spectrum is not 
scale invariant, which is explained in the text. 



action is ([L3|) . Thus it is found that 
3X 2 



^ + Hp) 



8X 4 



M u Ml 
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2M A {U - i) 4 

IhX 2 



l2Hip 3 
M 7 



12X 4 
M U M% 



2M 4 (t*-t) 4 



(25) 



(26) 



for Al(t* — t) 3> 1. In [26|, the results are applied to 
that of inflation, however, which are actually general for 
arbitrary evolution. Thus Q is given by 



Q = 



TX 



M 2 P {H 



M 7 Ml ' 



M 14 M 2 p F 

ip & 



M 2 Mp(U-t)\ 
(27) 

where Eqs.(QDl and are applied. Thus Q ~ |e| > 0, 
which is just required here, satisfies Eq.©. There is not 



5 



the ghost instability. Here, the importance of X 2 Uip is 
obvious, because if it disappears in (|15p. T is given by 



3X< 



1 



2M 4 (U - ty 



< 0, (28) 



Q > will hardly be obtained, which is consistent with 
Q = e < for this case. This indicate that it is X 2 Oip 
that alters the sign of Q, and leads Q ~ |e| > 0. The c 2 
is given by 

c 2 s = | ~ 1.4. (29) 

Thus > is constant, which is also just required. The 
sign of c 2 is determined by the signs of T and Q, both 
are positive. Here, obviously J- > is also required to 
assure c 2 > 0. Thus there are not the ghost and gradient 
instabilities, the effective theory is healthy. 

We plot the evolution of the amplitude of the curvature 
perturbation in Fig. 4, and the spectrum of perturbation 
in Fig. 5. We can see that the perturbation is initially 
not increasing, since it is inside the 1Z horizon. The in- 
crease begins until the perturbation mode leaves the 1Z 
horizon. The longer the wavelength of perturbation is, 
the earlier the perturbation leaves the 1Z horizon, the ear- 
lier it begins to increase. However, since the shorter the 
wavelength of perturbation is, the larger its initial ampli- 
tude is, all perturbation modes will eventually have same 
amplitude. 

There is a cutoff k cuto ff in Fig. 5, which is given by 



k, 



cutoff 



u 



inifr •> 



(30) 



where Wi n ifr is H freeze at initial time, and can be 
changed with the difference of the initial parameters in 
the numerical calculation. The spectrum is scale invari- 
ant for k > k cuto ff- However, for k < kcutoff > since 
the corresponding perturbation modes are outside the 7Z 
horizon all along, only are their amplitudes increasing but 
not the shape of the spectrum is not altered [4l| , [42| ■ 

The spectrum of 1Z is scale invariant. The amplitude 
of spectrum is given by Eq. p^|) 
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FIG. 6: The sketch of the evolution of the energy density p 
for different reheating courses discussed here. 



before this, the perturbation mode has leaved the Hubble 
horizon. 

Here, in certain sense, the reheating is alike with that 
for inflation. The preheating theory after inflation has 
been developed in 43] , 44] . In general, during the pre- 
heating phase after inflation the energy of inflaton will 
be rapidly released by the parametric resonance effects, 
due to the coupling of inflaton with other fields. Then 
this issue has been extensively studied, see (45j . (irij . (47| 
for reviews. 

We will apply the instant preheating mechanism [48| 
for given case here. We consider the straight coupling of 
tp with x particle as 



C ~ g 2 {ip - lfreh) 2 X 2 



(32) 



where g is the coupling constant. The effective mass of 
X particle is M 2 eff ~ g 2 (tp - ip re , 



. When the <p field 

arrives at the region around ip re h, M^eff ~ M xe ff, the 
adiabatic condition is broke, and the productions of x 
particles will be inevitable. This generally occurs in a 
region around ip re h, < ip re h/g, in which ip re h is the 
velocity of (p through ip re h- Thus the productions of \ 
particles is instantaneous, At re h ~ \j \J g(p r eh- 
The number density n x of x particle is 



1/2 

n 



M 



1/2 



(31) 



where A* ~ y/MMp is applied. Thus ~ 1U " re- 
quires M - 10- 10 c s M P . Thus M - 10 9 Gev for c s ~ I. 
The only adjusted parameter in this model is fixed by 
the observation. There is not other finetuning. 



I 

2^2 



f /2 €'el 

8tt 3 ! 



(33) 



where ilk is the occupation number of x particle. Thus 



Px = n x M x 



'Vrehi sinCe M xeff ~ 9{<P ~ Vreh) 



giprehAtreh. Thus the energy drained by the production 
of x particle is 



C. The reheating 



P<preh OK 



T ^ 

^reh ) 



(34) 



When the slowly expanding phase ends, the energy of 
Galileon field is required to be released into the radia- 
tion, and the universe reheats. Hereafter, the evolution 
of hot "big bang" cosmology begins. We can notice that 



where Eqs. ([T9|) and (|20l) arc applied, and p vre h is the 
energy density of tp around t re h- We assume tf ~ t re h 
for simplicity, i.e. the reheating occurs at the time when 
the slow expansion ends. Thus A4 2 Mp(t* — t re h) 4 ~ 1. 



This implies 



Px 

Pipreh 



g 2 M 2 
8tt 3 MJ, ' 



(35) 



We generally require M. <C 1 and g < 1. Thus 
p x /p V reh *C 1, which indicates that for such a single 
preheating, the energy of tp can hardly be released com- 
pletely, the universe is still dominated by p v , which will 
continue all along, since the energy density of tp is in- 
creasing with the expansion of universe while that of \ 
particle is decreasing. 

However, there might be M couplings, one of which is 
alike with (|3"2j) . We can find, after doing similar calcula- 
tions, that when 




FIG. 7: The figure of the effective potential for the exiting 
from the slow expansion. The black solid line is the motive 
trajectory of field in (ip, ip) space. 



TV > 



Mp 
g 2 M 2 



(36) 



the release of the energy of tp will be complete. The 
sketch of this reheating course is plotted in upper panel 
in Fig. 6. We assume that the x particle produced is 
rapidly transferred into the radiation. In this case, the 
reheating temperature T r is approximately determined 
by Pipreh ~ Tip . Thus we have 



•,10 



M 14 M 2 



1/4 



M^ 4 M 3 P /4 , 



(37) 



where M 2 Mp{t f — t re h) A ~ 1 is applied again. Thus if 



M ~ lCr lu M P , we have T r ~ 10 15 Gev. 

Here, TV 3> 1 is feasible, however, might be uncomfort- 
able. J\f 1 is required is because the energy of tp. has to 
be released completely for one time, or since the energy 
density of <p is increasing, the universe will dominated 
by tp all along. However, we also could consider another 
channel of the reheating, likes that in phantom inflation. 
The energy of ip is firstly shifted to the kinetic energy of 
a normal field, e.g.ip, and then the energy of ip is released 
by the instant preheating. The sketch of this reheating 
course is plotted in lower panel in Fig. 6. Here, the en- 
ergy of ip is not required to be released completely, since 
Pti ~ 1/a 6 is decreasing faster than that of the radiation, 
the universe will be dominated by that of the radiation 
early or late. 

We can implement it by considering a potential of tp, 
illustrated in Fig. 7. We require that it is only significant 
around or after |e| ~ 1, and is negligible |e| 3> 1. Then 
we introduce a waterfall field ip, coupled to tp. The effec- 
tive mass of ip is initially positive and becomes negative 
around |e| ~ 1. Thus ip will roll down along its potential. 
Thus almost all energy of tp will be shifted to p^, ~ ip . 
This energy will be expected to be released by the instant 
reheating. Thus there could be a suitable reheating after 
the slow expansion ends, after which the evolution of hot 
"big bang" cosmology begins. 



IV. DISCUSSION 

When initially e <C — 1 and is rapidly increasing, the 
universe is slowly expanding. The spectrum of curva- 
ture perturbation generated during such a phase of slow 
expansion can be scale invariant. This provides a mech- 
anism by which an alternative scenario of early universe 
can be imagined. Here, we show a model of such a 
scenario by applying an effective action of generalized 
Galileon. 

In principle, e < implies the ghost instability. 
However, in this model, because of the introduction of 
Galileon field, there is not the ghost instability, the per- 
turbation theory is healthy. In Refs.[l(|,[l9|, the phan- 
tom was applied for an implementing of slow expansion. 
In the calculations of perturbation, for consistence, |e| 
is used, though the initial value of perturbation is still 
pathologically defined. However, in the model given here, 
it can be found that actually Q ~ |e|. This in certain 
sense validates the argument and calculations used in 
[Tol| . fl9lj . i.e. the phantom field might be a simple sim- 
ulation of a full theory without the ghost below certain 
physical cutoff, which can give same results with that of 
a full theory, when the replacement of e with |e| is done. 

When e ~ — 1, the slow expansion ends. The exiting 
to a hot universe is only a simple reheating, since the 
universe expands all along. Thus there is not the prob- 
lem how the bouncing is implemented in bouncing cos- 
mologies [HI , i [H| • We have discussed possible imple- 
ments of reheating, and found that the available energy 
of Galileon field can completely released, the universe can 
reheat to a suitable temperature. Thus the model of the 
slow expansion given here might be a viable design of the 
early universe. 

The material compares of model with the observations 
is certainly interesting, which will place rigid constrains 
for the model. The results obtained will be expected to 
either improves or rules out this model. We will inves- 
tigate it elsewhere. However, it should be pointed that 
we only bring one of all possible implements of the slow 
expansion. In principle, there might be other effective ac- 
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tions of generalized Galileon, or modified gravity, which 
could give the same evolution of background. Thus for 
the slow expansion, it might be also significant to find 
alternative implements to the model given here, which 
will help to uplift the flexility of the slow expansion to 
the observations. 

Here, the scale factor is asymptotic to a constant value 
in infinite past, there is not singularity point. Thus in 
certain sense, the slow expansion scenario brings a solu- 
tion to the cosmological singularity problem. However, 
it also can be imagined that after the available energy 
of the field is released, it might be placed again in the 
bottom of its effective potential, and after the universe 
undergoes the radiation and matter periods, the field 
might dominate again and roll again with increasing en- 



ergy. This models an eternally expanding cyclic universe 
52], 53], 54], i.e. H oscillates periodically while a ex- 
pands all along. The implement of this cyclic universe 
might be interesting for refining with the model given 
here. 

Here, c s is constant is set. However, its change will 
obviously enlarge the space of solutions of the scale in- 
variance of curvature perturbation [55| . [561] . [57| . [58| . [59| . 
In certain sense all possibilities of the changes of a, Q 
and (? s might be interesting for further exploring. 
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